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A  DYNAMIC  MODEL  OF  A  REINFORCED 
THIN  PLATE  WITH  RIBS  OF  FINITE  WIDTH 


1.  INTRODUCTION 


Reinforced  plates  and  shells  have  a  variety  of  applications.  For  instance,  they  are  present 
in  the  design  of  ships,  undersea  vehicles,  and  aircraft  because  reinforcement  provides  increased 
structural  strength  with  minimal  increased  weight.  While  reinforcing  a  structure  will  allow  it  to 
survive  in  environments  in  which  an  unreinforced  structure  will  fail,  the  resulting  dynamic 
response  of  the  structure  changes  dramatically.  Early  work  in  the  area  of  reinforced  plates1'6 
generally  modeled  the  reinforcement  of  the  plate  as  a  line  stiffener  using  a  Dirac  delta  function  to 
mathematically  represent  the  effects.  Mace1'-  modeled  the  response  of  periodically  stiffened 
fluid-loaded  plates  to  harmonic  loading  and  to  line  and  point  force  loading.  Mace’s  work 
involved  a  transformation  into  the  wavenumber  domain  and  the  evaluation  of  a  contour  integral. 
Stepanishen  modeled  the  scattering  characteristics  of  a  plate  with  line  impedance 
discontinuities.  To  evaluate  the  scattered  pressure,  he  used  a  Fourier  integral  to  calculate  the 
relationship  between  the  plate  velocity  in  the  wavenumber  domain  and  the  spatial  domain. 

Eatwell  and  Butler4  determined  the  vibration  sound  radiation  from  a  fluid-loaded  plate  stiffened 
by  a  finite  number  of  beams.  T  heir  solution  involves  an  asymptotic  evaluation  of  the  pressure 
field  by  two  Fourier  integrals.  Cray^  determined  the  response  of  a  sectionally  aperiodic  plate  to  a 
line  force  in  the  wavenumber  domain.  Recently,  Hull6  derived  the  elastic  response  of  a  thick 
plate  system  to  harmonic  loading  with  stiffeners.  Some  work  exists  where  the  stiffeners  are 
modeled  with  finite  width.  Woolley  '  modeled  the  acoustic  scattering  from  a  plate  reinforced  by 
a  single  rib  and  by  a  finite  number  of  ribs.  In  these  papers,  Woolley  formulates  the  problem  in 
the  wavenumber  domain  and  then  solves  it  using  a  complicated  method  of  contour  integration. 
Woolley  allows  ribs  of  finite  width  in  his  model;  however,  he  states  that  another  author 
(Stepanishen)  has  obtained  a  different  result."  The  specific  problem  of  a  spatial  domain  response 


of  a  thin  plate  reinforced  by  ribs  of  finite  width  subjected  to  convective  loading  has  not  yet  been 
addressed. 

This  report  derives  an  analytical  model  of  a  thin  plate  of  infinite  spatial  extent  stiffened  by 
ribs  that  have  a  finite  width.  The  model  developed  here  is  different  from  previous  models  as  it  is 
based  on  differential  equation  theory  in  the  spatial  domain  only  and,  thus,  does  not  involve  a 
complicated  integral  to  be  evaluated.  The  governing  equation  is  a  flexural  plate  model  that  has 
an  external  load  and  is  reinforced  by  an  infinite  number  of  equally  spaced  ribs.  The  Heaviside 
step  function  is  used  to  load  the  rib  forces  onto  the  plate.  A  Fourier  series  then  replaces  the 
Heaviside  step  function,  and  it  is  shown  that,  with  this  substitution,  the  equation  decouples  using 
an  orthogonalization  procedure.  The  resulting  system  can  be  represented  by  an  infinite  set  of 
algebraic  equations.  These  equations  are  truncated,  and  a  solution  to  the  plate  displacement  is 
found.  An  example  problem  is  formulated  and  the  results  are  compared  to  finite  element  theory 
to  ensure  that  the  proper  analytical  solution  has  been  obtained.  The  results  are  discussed. 
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2.  SYSTEM  MODEL  AND  ANALYTICAL  SOLUTION 


The  equation  governing  the  motion  of  a  plate  with  an  applied  external  foree  and  an  infinite 
set  of  finite  width  stiffeners  ean  be  derived  through  a  force  balanee  along  the  length  of  the  plate. 
This  differential  equation  written  in  the  spatial-time  domain  is 


-34k<jc./)  ,d2w(x,t)  x  ....  ..  /  .u/  .  fx 

D-  —  +  ph - - —  =  -f(xJ) /  w(x,t)\\(x- nL)+ —  /  /)H(.v  —  b  —  nL) , 

A  A 


d.\- 


dr 


n=-*j 


(i) 


where  w(x,l)  is  the  transverse  displacement  of  the  plate  in  the  v-direction  (m),  x  is  the  spatial 
location  on  the  plate  (m),  t  is  time  (s),  p  is  the  density  of  the  plate  (kg  m  '),  h  is  the  height  of  the 
plate  (m),/(.\\/)  is  the  applied  external  load  on  the  plate  (N  nf‘),  K  is  the  stiffness  of  eaeh  rib  per 
unit  length  (N  nf2),  L  is  the  distance  between  the  left  edges  of  adjacent  ribs  (m),  b  is  the  width  of 
each  rib  (m),  H  is  the  Heaviside  step  function  (dimensionless),  and  D  is  the  flexural  rigidity  of 
the  plate,  which  is  given  by 


Eh3 

12(1  -v1 )  ' 


(2) 


where  v  is  Poisson’s  ratio  of  the  plate  (dimensionless)  and  E  is  Young's  modulus  of  the  plate 
(N  nf2).  In  equation  (1 ),  the  first  term  with  the  Heaviside  function  corresponds  to  the  left  edge 
of  the  ribs  and  the  second  term  corresponds  to  the  right  edge  of  the  ribs.  A  diagram  of  the 
system  subjected  to  a  continuous  spatial  force  is  shown  in  figure  1 .  The  foreing  function  and  the 
response  are  harmonie  in  time  (i.e.,  J (xj)  =  F(x)exp(-\(ot)  and  w(x,l)  =  W’(.Y)exp(-iru/)); 
thus,  equation  ( 1 )  can  be  written  in  the  spatial-frequeney  domain  with  an  applied  foree  at  a 
definite  wavenumber  as 


d4HYvl  K  "~+y  K "~+ 00 

D - -  phw2W(x)  =  —Fq  exp(iA.v)  -  —  (a  )H(.y  -  nL)  +  —  ^lT(.v)H(a  -/i-/;Z.)  , 


The  Heaviside  step  function,  H,  also  called  the  unit  step  function,  is  a  discontinuous  function  whose  value  is 
zero  for  negative  argument  and  one  for  positive  argument.  It  seldom  matters  what  value  is  used  for  H(0).  since  H  is 
mostly  used  as  a  distribution.  The  function  is  used  in  the  mathematics  of  control  theory  and  signal  processing  to 
represent  a  signal  that  switches  on  at  a  specified  time  and  stays  switched  on  indefinitely. 
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Figure  1.  Schematic  of  Reinforced  Plate  with  an  Applied  Spatial  Force 

where  F{)  is  the  magnitude  of  the  applied  force  (N  m'~)  and  A  is  the  wavenumber  of  the  applied 
load  (rad  m  ').  Because  of  the  spatial  periodicity  of  the  system,  the  magnitude  of  the  response 
can  be  expressed  as9 


where 


/7?=+00 

W(x)=  Qxp(ikmx)  , 

00 


(4) 


km  =  k  + 


l/rm 


(5) 


so  that  A0  s  k  ,  and  the  Wm's  are  the  unknown  coefficients  whose  solutions  are  sought. 
Substituting  equation  (4)  into  equation  (3)  and  evaluating  the  derivatives  yields 


7/?=+00 


X  ( Dkm  “  ph(°2  W,„  exp(iAnix) 

/;=+x 


;;;=— oo 


K 


=  ~Fq  exp(ifcv)  +  —  Y 

b 


jj;=+oc 


exp(iA-„rv) 


[H(x-h-nL)-H(x-nL)]  , 


which  is  an  algebraic  summation  problem  free  of  differentials. 


(6) 
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The  presence  of  the  Heaviside  step  functions  in  the  third  term  complicates  the  form  of 
equation  (6).  However,  if  the  Heaviside  functions  are  expressed  as  a  Fourier  series,  not  only  are 
the  difficulties  associated  with  these  discontinuous  functions  eliminated,  but  also  each  term  of 
the  equation  becomes  an  exponential  function  of.v.  The  Fourier  series  of  the  Heaviside  functions 
is  written  as 


77 -+-TC  II— -hOC 

V[H(.y  -  nL)  -  H(.y —  b  —  nL)]  =  YjCn  exp(i2/r/7.v/  L) , 

/;—oc  n= — oo 


(7) 


where 


exp(-i 2xnbl  L) 
\2nn 


h 

L 


n  it  0 


//  =  0 


(X) 


Inserting  this  Fourier  series  identity  into  the  third  term  of  equation  (6)  and  extracting  functions 
independent  of  n  and  m  from  the  double  summation  yields  the  term 


is  77-+00 

- 1 
b  ~ 


/77=+X) 


eXP('kMX) 


[H(x  -b-  nL)  -  H(.v  -  nL)\ 


77=+0O 


=  -  — cxp(i^Y)  Y 


777-+0O 


Y^'m  CXp(i2/T»?.Y/L) 


c„  exp(i2/r«x/L) . 


(9) 


Furthermore,  the  double  summation  in  equation  (9)  can  be  rewritten  as 


77=  +X 


exp(ifor)  Y 


K 


77= — QT- 
77=  TX  777=  fX 


/77=+X 


Y^'w  exp(i2^WA :/L) 


c„  c.xp(\2iznx/L) 


=  ~ ~  Y  Z  Wn  Cm  v  CXP('A  „,*)  ■ 


7/=-X  777= — X 


(10) 


The  proof  of  the  double  summation  identity  in  equation  ( 10)  is  presented  in  the  appendix. 
Inserting  the  second  part  of  equation  (10)  into  equation  (6)  yields 
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777=+X 


K 

X ( Dkm  -  P]uo~  VV»>  exP('A'mx)  =  ~Fo  exp(ifcx-)  -  —  X  X  W»  c>«-»  exp(iA„rv)  ,  (11) 

7/7— — CC  }  U--rr  777“— 00 


77  — ~HX  777=+X 


whieh  is  the  algebraic  equation  that  models  the  stiffened  plate  subjected  to  an  applied  foree. 

The  solution  to  the  unknown  W,„  coefficients  is  now'  found  by  an  orthogonal  expansion  of 
equation  (11).  Specifically,  equation  (1 1)  is  multiplied  by  the  exponential  exp(-iA^A')  and  the 

resulting  expression  is  integrated  over  [0,L].  Beeause  the  exponential  functions  exp(-ik px)  and 

exp(iA/77.v)  are  orthogonal  on  this  interval  when  m  ^  p ,  equation  (11)  decouples  into  an  infinite 
number  of  individual  /^-indexed  equations,  each  one  expressed  as 


( Dp, -pi, a,2  W ,, 


/?=+00 
Xi  />-» 

77— —X 


W„ 


[~F{)  P  =  0 

I  0  0 


(12) 


Writing  out  all  the  equations  from  equation  (12)  for  -oo<  p  <ooand  plaeing  them  into  matrix 
form  yields 


A_  |  0  0 
0  ‘  A0  0 
0  0/1, 


• 

w_, 

W0 

K 

)  +  — 

b 

IV, 

c0  c-l  c-2 

c\  c0  c-\ 

c2  c’l  c0 


: 

* 

W-, 

0 

W0 

ll 

—  Fo  > 

W\ 

0 

: 

t..  *  i 

(13) 


where 

A/t  =  Dkp  -  phco 2  . 


(14) 
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The  first  term  of  equation  (13)  represents  the  dynamics  of  the  plate,  the  second  term  represents 
the  dynamics  of  the  ribs,  and  the  third  term  represents  the  external  load  on  the  structure.  Note 
that,  mathematically,  the  effect  of  the  ribs  is  the  constant  (K/b)  multiplied  by  a  matrix  that 
contains  a  pemiutation  of  the  Fourier  coefficients  of  the  Heaviside  step  function.  The  unknown 
coefficients  Wm  in  equation  (13)  can  be  solved  by  truncating  the  matrices  and  vectors  to  a  finite 
number  of  terms  and  analytically  evaluating  the  equation 


I 

• 

-1 

A_  i 

0 

0 

c0 

c_\ 

£•_-> 

0 

K 

"o 

►  = 

...  0 

4) 
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4-  — 
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£■_,  ••• 
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~Fo  > 

0 

0 

4 

1 1 

C’2 

c\ 

c0 

0 

L  *  J 

• 

(15) 


3.  NUMERICAL  EXAMPLE 


A  numerical  example  is  now  formulated  and  discussed.  The  problem  consists  of  a  thin 
plate  that  has  a  Young’s  modulus  of  3  x  10*  N  m ",  density  of  1200  kg  m  \  Poisson’s  ratio  of 
0.45,  and  a  height  of  0.01  m.  Each  rib  has  a  stiffness  per  unit  length  of  1  x  10s  N  m  2  and  a  width 
of  0.2  m.  The  center-to-center  spacing  of  adjacent  ribs  is  l  m.  The  applied  load  has  a  magnitude 
of  1  N  m  2  and  a  wavenumber  of  0  rad  m  '.  These  parameters  are  input  into  equation  (15)  to 
determine  the  coefficients  that  are,  in  turn,  entered  into  equation  (4)  to  calculate  the  displacement 
of  the  plate.  Figure  2  is  a  plot  of  the  magnitude  of  plate  displacement  in  the  decibel  scale  versus 
the  spatial  location  at  a  frequency  of  50  Hz.  The  solid  line  is  the  solution  derived  in  this  report, 
while  the  dots  are  a  verification  solution  determined  using  the  finite  element  method.  The 
dashed  line  on  the  plot  indicates  the  spatial  position  of  the  right  edge  of  the  rib  that  corresponds 
to  n  -  0  in  equation  ( 1 ).  This  specific  problem  converged  using  3 1  terms  (  -  1 5  <  p  <  1 5 ), 
although  it  is  noted  that  convergence  is  a  function  of  the  specific  parameters  used  in  the  model. 

Several  features  in  figure  2  are  now  noted.  First,  the  platc/rib  system  is  softer  between  the 


ribs  and,  thus,  has  a  greater  displacement  at  these  spatial  locations.  As  a  result,  if  the  stiffness  of 
the  ribs  approaches  0,  the  solution  approaches  that  of  an  unreinforced  plate.  Second,  even 
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though  the  ribs  are  fairly  stiff  in  this  example,  they  still  exhibit  dynamic  behavior  aeross  their 
spatial  extent.  As  evidence,  examine  figure  2  in  the  area  from  x  =  0  to  x  =  b  =  0.2.  Each  rib  has 
almost  40  dB  of  motion  in  this  relatively  small  spaee,  although  the  displacement  is  significantly 
lower  than  the  plate  displacement  between  the  ribs.  Third,  although  not  explicitly  shown  in 
figure  2,  the  solution  is  periodic  on  „v  e  [0,  L ] .  Seetional  aperidoeity  can  be  added  to  the  model 

using  previously  developed  techniques.' 


Figure  2.  Magnitude  of  Plate  Displacement  Versus  Spatial  Location  Showing  the  Analytical 
Solution  ( — )  and  the  Finite  Element  Solution  (•) 

(The  dashed  line  is  the  location  of  the  right  edge  of  the  n  =  0  rib.) 


8 


4.  CONCLUSIONS 


A  model  of  a  reinforced  thin  plate  with  ribs  of  finite  thickness  was  derived  and  the  solution 
was  obtained  using  an  orthogonalization  procedure.  It  was  shown  that  the  effect  of  the  ribs  is, 
mathematically,  a  constant  multiplied  by  a  matrix  that  contains  a  permutation  of  the  Fourier 
coefficients  of  the  Heaviside  step  function.  A  numerical  example  was  presented  and  compared 
to  a  finite  element  analysis.  This  new  analytical  solution  was  in  almost  exact  agreement  with  the 
finite  element  solution.  The  results  of  the  model  and  the  associated  dynamics  were  briefly 
discussed. 
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APPENDIX 

DOUBLE  SUMMATION  IDENTITY 


The  proof  of  the  double  summation  identity  given  in  equation  ( 10)  is  presented  in  this 
appendix.  The  identity  is  critical  in  transforming  the  problem  into  a  tractable  algebraic  equation. 
The  first  term  in  equation  (10)  without  the  (K/b)  coefficient  is 
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777=+0O 
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To  simplify  the  notation  of  the  following  equations,  the  formula 

=  exp(i2/r/77.Y/L)  (A-2) 


is  used.  This  term  has  the  property 
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Expanding  the  m  indexed  series  results  in 
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Next,  expanding  the  n  indexed  series  yields 
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and  multiplying  through  gives 
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Regrouping  equation  (A-6)  on  specific  values  of  ff/7  yields  the  equation 
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which  can  be  rewritten  as 
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